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We show that the semiclassical dynamics of an electronic droplet confined in the plane in a 
quantizing inhomogeneous magnetic field in the regime when the electrostatic interaction is negligible 
is similar to viscous (Safi'man- Taylor) fingering on the interface between two fluids with different 
viscosities confined in a Hele-Shaw cell. Both phenomena are described by the same equations with 
scales differing by a factor of up to 10~^. We also report the quasiclassical wave function of the 
droplet in an inhomogeneous magnetic field. 
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An important class of pattern formation on moving 
fronts occurs when diffusion, rather than convection, 
dominates the transport. In these cases a front is driven 
with the normal velocity proportional to the gradient of 
a harmonic field - a mechanism known as D'Arcy's law 
or Laplacian growth (for a review see, e.g., Viscous 
or Saffman- Taylor fingering is one of the most studied in- 
stabilities of this type. It occurs at the interface between 
two incompressible fluids with different viscosities when 
a less viscous fluid is injected into a more viscous one 
in a 2D geometry (typically, the fluids are confined in a 
Hele-Shaw cell - a thin gap between two parallel plates) 
or in porous media The interface forms a pattern 
of growing fingers whose shape becomes complex at high 
flow rates (Fig. 1). 

Instabilities on diffusion driven fronts occur in differ- 
ent settings ranged from geological to molecular scales. 
In this letter we show that similar phenomena may take 
place in semiconductor nanostructurcs. A growth of an 
electronic droplet in a Quantum Hall (QH) regime, while 
increasing the number of particles in the droplet, is sim- 
ilar to Laplacian growth. 

We first recall the Saffman- Taylor (ST) problem. In a 
thin cell, the local velocity of a viscous fluid is propor- 
tional to the gradient of pressure: v — — Vp (the D'Arcy 
law). Incompressibility implies that the pressure p{z) is 
a harmonic function of z = x + iy with a sink at infinity: 



V2p(z) = 0, piz) 



\og\z\, |z|^oo (1) 



If the difference between viscosities is large, the pressure 
is constant in the less viscous fluid and, if the surface 



•Electronic address: igamOamadeus. fiz.huji.ac.il 



tAlso at Landau Institute of Theoretical Physics 

tAlso at ITEP, Bol. Cheremushkinskaya str. 25, 117259 Moscow, 

Russia 



tension is ignored, it is also constant (set to zero) on the 
interface C. Thus on the interface 



p{z) = 0, 



-dnP{z) 



zeC. 



(2) 



At constant flow rate, the area of the less viscous fluid 
grows linearly with time t. We set the flow rate such that 
the area equals irt. The other parameters, viscosity and 
the width of the cell, are chosen to set pressure equal to 
the velocity potential. For recent advances in the studies 
of flngering in channel and radial geometries see ^. 

It is convenient to rewrite Eqs. (|l|,^) in terms of the 
conformal map, w{z,t), of the domain occupied by the 
more viscous fluid to the exterior of the unit disc |it;| > 1 
in such a manner that the source at z = oo is mapped 
to infinity. In terms of the conformal map the pressure 
is p ^ —^\og\w{z,t)\ and the complex velocity in the 



viscous fluid is v{z) 



IV y 



\dz\-Ogw{z). On the 



interface, it is proportional to the harmonic measure: 



Vn{z,t) = ^\w'{z,t)\, 



z€C. 



(3) 



The complex velocity is conveniently written using the 
Schwarz function, S{z), 

dz\ogw{z) = dtS{z). 

The Schwarz function is analytic in a domain containing 
the contour C such that S{z) — z on C 

In the idealized settings (0J|), the Saffman- Taylor 
problem confronts an obstacle. As a result of the scale 
invariance, some fingers develop cusp-like singularities 
within a finite time |^. A modification of the growth 
law which introduces a mechanism curbing the curvature 
of the interface at a microscale is necessary. 

The known cut-off mechanisms (e.g., surface tension, 
lattice regularization, etc.) destroy the mathematical 
structure of the idealized problem and make it difficult 
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FIG. 1: A schematic illustration of a less viscous fluid domain 
or an electronic droplet in a strong magnetic field. Electronic 
droplet is stratified by semiclassical orbits. The domain of a 
more viscous fiuid (exterior of the droplet) mapped confor- 
mally into the exterior of a unit circle. 

for analytical analysis. It is believed, however, that once 
a steady fractal pattern has been developed, its fractal 
character does not depend on the mechanism stabilizing 
the singularities. 

The purpose of this letter is twofold. One is to show 
that a natural mechanism which introduces a scale but 
captures the physical (and the mathematical [Q) struc- 
tures of the problem is "quantization". Quantization im- 
plies that a change of the area of the domain is quan- 
tized. This ensures that no physical quantity can have 
features on a scale less than a "Planck constant" which 
cuts-off singularities. The Saffman- Taylor problem (jljj^) 
then arises as a semiclassical limit, when the scale intro- 
duced by "Planck's constant" tends to zero. 

Our second goal is to show that quantized Saffman- 
Taylor problem describes an electronic droplet in a spe- 
cial QH-regime where electrostatic forces are weak. This 
correspondence suggests that the edge of the QH-droplet 
may develop unstable features similar to the fingers ob- 
served in Hele-Shaw cell (a fingering pattern is illustrated 
in Fig. 1). A spatial resolution of nanoscale structures in 
GaAs by recently developed scanning techniques might 
be apt for a search for such features H. 

We will study a shape of a large electronic droplet 
on the fully occupied lowest Landau level of a quantiz- 
ing magnetic field. The magnetic field is assumed to be 
nonuniform in the area away from the droplet. We show 
that Aharonov-Bohm forces, associated with the nonuni- 
form part of the magnetic field shape the edge of the 
droplet in a manner similar to a fingering interface driven 
by a Laplacian field. 

In order to present our argument we shall neglect the 
interactions among the electrons and assume that the ex- 
ternal electrostatic potential is zero. Under these condi- 
tions we will show that the semiclassical dynamics of the 
QH droplet is governed by the same equations of viscous 
fingering scaled to a nanometer scale. By the semiclassi- 
cal limit we mean a large number of electrons iV, small 
magnetic length but a finite area of the droplet: 

TV ~> oo, £= yJhc/eBo ^0, t = 21^ N = finite. 

The droplet grows when its quantized area changes by 
a quantum, 7ri — > 7ri -f 27rf^, as an electron is added 
{N ^ TV 4" 1) or by changing the magnetic length. 

Prior to computations, we would like to give an insight 
why electronic droplet in magnetic field may behave sim- 



ilarly to viscous fingering. Equations (|i|, H) imply that 
the harmonic moments of the viscous fluid domain, 

tk = -\ / z-'^d'z, fc = l, 2,... 

'^'^ Jvisc. fluid 

do not change in time They are initial data 

of evolution. Indeed, ^t^ = §^ z^'' dnp{z)\dz\ = 
<f^ z~''dz \ogw{z)dz = since dz log w(z) is analytic 
in the viscous domain. Conservation of the harmonic mo- 
ments is an equivalent formulation of D'Arcy's law (^. 

As is shown below, the harmonic moments of the QH- 
droplet do not depend on the number of particles in the 
droplet. The moments feature the inhomogeneity of the 
magnetic field. A proof of this assertion would yield 
D'Arcy's law for the QH droplet. 

Our proof consists of two steps. First we construct the 
wave functions of the QH-states in a weakly nonuniform 
magnetic field and introduce the r-function - a generat- 
ing function of the matrix elements. Then we determine 
semiclassical states - the orbits (supports of the wave 
functions) , and study an evolution of the orbit as the area 
of the droplet is increased. We will consider the droplets 
with smooth shape, leaving analysis of the finger-tip sin- 
gularities for further publications. The size of the letter 
forces us to rely on some formulas on the evolution of 
conformal maps proved in pO| | . 

Let us first recall the concept of a QH-droplet (see e.g., 
|l^). Consider spin polarized electrons on a plane 
in the lowest level of a quantizing nonuniform magnetic 
field, directed perpendicular to the plane, B[x,y) > 0: 

The lowest level of the Pauli Hamiltonian is degenerate 
even for a nonuniform field. The degeneracy equals the 
integer part of the total magnetic flux $ = J dxdyB in 
units of flux quanta, $0 = 2Trh (we set e = c — 1). The 
orthonormalized degenerate states, written in transversal 
gauge, are ipnixju) — P„(z)exp-^, where the potential 
W{x,y) obeys the equation B — —\/^W, and Pn{z) are 
holomorphic polynomials of a degree n, which does not 
exceed the degeneracy of the level . 

We will consider the following arrangement: A strong 
uniform magnetic field -Bo > is situated in a large disk 
of radius Rq; The disk is surrounded by a large annulus 
Rq < \z\ < Ri with a magnetic field i?i < directed 
opposite to Bq, such that the total magnetic flux <& of 
the disk, \z\ < Ri, is N^q. The magnetic field outside 
the disk \z\ < Ri vanishes. The disk is connected through 
a tunneling barrier to a large capacitor that maintain a 
uniform small positive chemical potential slightly above 
the zero energy of the lowest Landau level. 

In this arrangement a circular droplet of N electrons 
is trapped at the center of the disk \z\ < Rq. We choose 
the magnetic field Bi such that its radius 1\J2N is much 
smaller than the radius of the disk i?o- 

Next we assume that a weakly nonuniform magnetic 
field 8B is placed inside the disk \z\ < Rq but well away 
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from the droplet. The nonuniform magnetic field does 
not change the total flux J SBdxdy = 0. The droplet 
grows when Bi is adiabatically increased, keeping Bq and 
6B fixed. Then the degeneracy of the Landau level and, 
consequently, the size of the droplet increase. 

The wave function of the droplet is the Jastrow func- 
tion 



1 



, zn) = 



VNWn 



A(z)e 



(4) 



where 



A(z) = Y\_ ~ ^rn) = det (P„(Zm)) 



n,m<N 



n<m<N 



is the Van der Monde determinant, and the normalization 
factor - the r-function- 



TN 



1 

TV! 



(5) 



depends on magnetic field. 

In terms of the orthonormal one-particle states, '0„(z), 
the density pn{z) = N J [^(z,^!, . . . , ^N-i)\'^d'^^ reads 

PN+i{z) - Pn{z) = \'4>n{z)\^. 

This formula defines a growth process. One-particle state 
gives the probability of adding an extra particle to a point 
z of the droplet. 

The integral in (|^) converges in the area of the droplet, 
where magnetic field is uniform. In this area the annu- 
lus, i?o < 1^1 < Ri, does not contribute to the potential 
W , and the potential V of the nommiform part of the 
magnetic field 5B — —V'^V is a harmonic function 

V{z) = i?e ^ = -i- /" 5B{z)z-''d^z. 

k>l '' 

The parameters tk are, now, the harmonic moments of 
the deformed part of the magnetic field. Summing up, 
we have W{z) = -\z\^l^e ^ V[z). 

The one-particle state ipNiz) can be expressed in terms 
of the T-function (further below we set 2£'^ — K) 



IpNiz) 



.1(1 



^V{z)) 



yJ{N + 1)tn+itn 



.^e-''^(^)r^., (6) 



where D{z) — X]fc>i ^irW' ^ generator of the defor- 
mations of the magnetic field. 

The next step is a semiclassical analysis of the wave 
function, ^n{z). In this limit the droplet has a sharp 
boundary (we assume that the boundary is smooth) , and 
the semiclassical states (with large N) are localized at 
orbits of width of the order of t. Since the magnetic field 
in the area of the droplet is assumed to be uniform, the 
density is also uniform inside the droplet. The area of 
the droplet is quantized as = 2ttNP . 



As ft — > with t = hN fixed, the integral (||) is sat- 
urated by its saddle point. At the saddle point, the ^„ 
are uniformly distributed in a domain determined by the 
area ttTlN and the harmonic moments tk jl^ . We will see 
that this domain determine the shape of the droplet. At 
the saddle point the r-function (||) tends to the classical 
T-function F = lim^j^o log T^r , studied in [nO[ : 



F{t;tu---) = 
On expanding 





1 1 


/ / log 




/ J droplet 


z z' 



d'zd'z'. 



and 



A;<1 



A{z,z) 



in h, and defining 
1 



t log z 



dt + D{z) F, 



one can write the amplitude of a semiclassical state, in 
terms of variations of the classical r-function: 



IV'iv(^)|' 



,-i(9?-2fle_D^(z))F 



^-iA(z,z)^ i = ftiV. (7) 



Properties of the functions n(z) and D^{z)F clarify the 
meaning of a semiclassical state. We list them below. 
The function A{z) is the imaginary part of the semiclas- 
sical action. It was proved in |l^ that A 

- reaches its minimum on a closed contour (an orbit) 
which bounds the domain with harmonic moments 
tk and the area 7rt, and everywhere on the orbit the 
action vanishes: 



A{z) = 0, Ren{z) = ^\z\ 



z G orbit. 



The first derivative of the action normal to this 
contour also vanishes on the orbit: 



d,A{z) = 0, d,n{z) ^z, ze orbit. 



(8) 



- If the area of the contour increases, while tk remain 
fixed, the action changes as 

dtA = -~dtRen{z,t) = - log\w{z,t)\. 

Eq.(||) implies that S{z) = dzfl{z) is the Schwarz func- 
tion of the contour. 

The the action A remains positive in the vicinity of 
the contour and everywhere in the exterior domain. Its 
second variation normal to the contour reads 

Aiz + 6nZ) = \SnZ\'^ + ... , zeC 

where SnZ is a normal deviation from a point z on the 
contour (we use the relation |5„zp = — S"(z)((5„z)^). 
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A semiclassical state is localized at the minimum of 
A{z), where the amplitude has a sharp maximum. A 
contour where the action vanishes is the orbit. We con- 
clude that all orbits have the same harmonic moments tk 
and differ by the area. This implies the D'Arcy law. 

The second variation of the classical r-function was 
also computed in Ref. | |lO| |: 

Merging the points z and z' , we find that the value of the 
first factor in is just the harmonic measure |w'(z)|. 

Summing up, the amplitude of the semiclassical wave 
function in a weakly inhomogeneous magnetic field is 

|^^(,)|2^M£le-|^(^,^-), t = Nh. 

TT\/2'Kh 

The factor in front of the exponent ensures the correct 
normalization: ji/jArpds = d0(27r'^/i)^-^/^e^R l'^"^' (here 
ds is an element of the length of the orbit, and dcfi is an 
element of the angle on the unit circle). The classical 
probability distribution is then IV'Arp — j^lw' (z)\Sc{z), 
where the 5-function is localized on the orbit. 

When a new particle is added to the droplet, the edge 
advances. The velocity Vn of the advance is defined as 
dtp{z) — Vndnp{z). The normal derivative of the den- 
sity is d„p — -^Sc{z), and a change of the density is 
hdtp{z) ~ pN+i - PN = IV'A'P - ^\w'iz)\Sc{z). This 
prompts D'Arcy's law (||). 

It is interesting to compare dimensionless viscosity of 
liquids used in viscous fingering experiments p^ , and a 



"viscous" effect of quantum interference of at the first 
Landau level. The parameter of the dimension of length 
controlling viscous fingers in fluids - an effective capil- 
lary number - is ^j^f, where q is the flow rate, b is 
the thickness of the cell, t] is the viscosity and a is the 
surface tension. In recent experiments [ |l6| , this number 
stays in high hundred of nanometers, but can be easily 
decreased by increasing the flow rate. This length is to be 
compared with the magnetic length £. At magnetic field 
about 2T it is about 50 nm. At these numbers, the cap- 
illary effects are less important in electronic liquids than 
in classical ffuids, but remarkably stay in the same range. 
Under conditions reducing electrostatic effects, fingering 
in electronic droplets may be even more dramatic than in 
viscous fiuids. Semiconductor devices imitating a chan- 
nel geometry of the original Saffman- Taylor experiment 
may facilitate fingering instability. 
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